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Probing the sheath electric field with a crystal 
lattice by using thermophoresis in dusty plasma 

Victor Land, Bernard Smith, Lorin Matthews, Truell Hyde Fellow, IEEE 
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Abstract — A two-dimensional dust crystal levitated in the 
sheath of a modified Gaseous Electronics Conference (GEC) 
reference cell is manipulated by heating or cooling the lower 
electrode. The dust charge is obtained by measuring global 
characteristics of the levitated crystal obtained from top-view 
pictures. From the force balance, the electric field in the sheath 
is reconstructed. From the Bohm criterion, we conclude that the 
dust crystal is levitated mainly above and just below the classical 
Bohm point. 

Index Terms — Dust crystal, dusty plasma, modified GEC cell, 
sheath electric field, thermophoresis. 



I. Introduction 

DUST particles present in laboratory plasma collect elec- 
trons and ions and become negatively charged. The 
charge on a particle is estimated to be between 1 and 5 
electron charges per nanometer radius [1], depending on 
the background pressure 0, 0. Therefore, micrometer- sized 
particles used in dusty plasma experiments can carry a negative 
charge of thousands of electron charges. Because of this, the 
dust particles are levitated against gravity in the strong elec- 
tric fields present in the sheath. The electrostatic interaction 
between the particles is shielded by the plasma, but depending 
on the discharge parameters, strongly coupled systems can be 
formed Q] . By applying a radial confinement, these so called 
dust crystals can be used to study many solid state physics 
phenomena on a scale accessible to ordinary optical diagnostic 
techniques. 

To understand the behavior of these systems, knowledge 
about the screening length and the dust charge is of paramount 
importance. Furthermore, once the charge on the dust particles 
is known, the vertical force balance between gravity and the 
electric field can provide insight into the structure of the 
electric field in the sheath. This is important for both dusty 
plasma experiments and many plasma enhanced industrial 
processes, where ions are accelerated by the sheath electric 
field from the plasma bulk to a processing substrate l5l . 

Different techniques using dust particles suspended in the 
sheath of a plasma as probes to measure properties of the 
plasma have been developed in the past. Some techniques 
depend on measuring the equilibrium levitation height of a few 
dust particles, while using passive emission spectroscopy, or 
probe measurements to obtain additional parameters 151. 171. 
Other techniques use perturbations of the plasma by adding 
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a low-frequency sinusoidal signal to the lower electrode, or 
sudden changes in the DC bias (8], or perturb the dust particles 
by using lasers to excite dust oscillations __|. 

Despite being successful, these methods do have several 
difficulties; often only a few dust particles are used at a time, 
limiting the spatial range of measurements at a given time, 
optical measurements in the sheath are technically challenging 
in the small volume close to the highly reflecting electrode 
ifTOlL while optical access can be hard to obtain altogether 
ifTTl . Probe measurements are usually too perturbative to dusty 
plasma to be reliable, while induced perturbations to the 
plasma inevitably cause changes to the dust parameters, such 
as the charge. 

By manipulating a two-dimensional dust crystal as a whole 
with thermophoresis by heating and cooling the lower elec- 
trode 1 12 ] and by measuring global observable quantities from 
top-view pictures, it is possible to reconstruct the electric field 
in the sheath by using an analytical model fl3l . Here, we 
present the experimental reconstruction of the sheath electric 
field above the lower powered electrode in a modified GEC 
cell fl4l . One of our conclusions is that the dust is mostly 
levitated on the plasma side of the classical Bohm point, 
consistent with previous findings. 

II. Experimental setup 

A standard 13.56 MHz radio-frequency (RF) discharge in 
argon was created in a GEC cell, which is modified for dusty 
plasma experiments. Cover plates with cylindrically shaped 
cutouts milled in them can be attached on top of the powered 
bottom electrode. These cutouts provide an electric field that 
confines the dust particles radially. The upper electrode is 
a hollow grounded cylinder, allowing optical access to the 
dust crystal from above. Above this electrode, dust containers 
are installed, which release dust particles when tapped. The 
particles then fall down through the hollow electrode into the 
plasma, until they reach their equilibrium levitation height. 

A camera mounted above the hollow upper electrode is used 
to provide top-view pictures of the horizontal dust crystal 
structure. A second camera installed in front of a side-port 
takes side-view pictures of the vertical structure of the dust 
crystal. Two red diode-lasers with wavelengths of 640 nm and 
686 nm are equipped with cylindrical lenses that shape the 
laser beams into thin laser sheaths. One sheath is oriented 
vertically, the other horizontally. The light scattered by the 
dust particles is captured by either the top-view or side-view 
camera. Filters can be used with the cameras to only allow 
light at the laser wavelength to be observed, rejecting the 
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plasma glow. Typically, only one camera-laser combination 
is used at a time. The cameras and lasers can be moved in 
three directions; however, the focus of the cameras is fixed, 
so that the resolution does not change between pictures. The 
resolution of the top-camera and side-camera were measured 
to be 26 and 14 fim per pixel, respectively. 

The lower electrode can be cooled or heated by running 
liquid through a circular channel in the electrode. For heat- 
ing and cooling, we used a MGW Lauda Brinkman RM3T 
chiller/heater, connected to a Ranco ETC-1 11000-000 ther- 
mostat. The temperature was measured with thermoresistors 
attached near the thermostat, as well as on the bottom of the 
electrode. Due to the hysteresis of the thermostat/chiller-heater 
system, the temperature fluctuated approximately 2° C above 
and below the set temperature. A rapid rise to temperatures 
above the set temperature was followed by a slow decrease 
to a point 2 degrees below the set-point. By measuring the 
electrode-temperature variation, we determined an interval of 
at least 60 seconds where the temperature change was slow and 
the temperature could be considered constant. This was long 
enough to take the necessary top- and side-view pictures, since 
the cameras operated at 60 and 30 Hz, respectively. Figure [T] 
shows the change in temperature of the lower electrode and 
the ambient temperature. 

To determine the electric field profile, the temperature of 
the lower electrode was varied from 0° C to 60° C in steps of 
10° C. Since the electrode and the cover-plate are both made 
of metal and are heated relatively homogeneously, we expect 
temperature gradients along the surface of the electrode to 
be small. We are therefore not concerned with thermal creep 
effects recently reported to play an important role in dusty 
plasma contained within glass walls lH5l . 
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Fig. 1. The room-temperature and the temperature of the powered electrode 
measured with thermoresistors. The room-temperature remains close to 23° C. 
The temperature on the electrode was set at 20° C. It can be seen that 
the heater/thermostat rapidly heats the system, but then overshoots by a few 
degrees. After reaching the maximum temperature, the electrode slowly cools 
down, giving a relatively long time-window where the temperature may be 
considered constant. 



For the experiments presented here, 8.89 /im diameter 
melamine-formaldehyde (MF) particles were used and the 
discharge was run at 2 W of input power and 200 mTorr 
background pressure. We used a 1-inch (25.4 mm) diameter 
cutout, being roughly 0.5 mm deep. 



III. Observables 

We assume that the radial confinement potential provided 
by the cutout is parabolic, so that we can write the potential 

as 



(r, z) = (j)(z — h(r)), h(r) = cr 2 



(1) 



which can be visualized as vertically stacked parabolically 
shaped equipotential surfaces with the central axis of the 
cutout (with r = 0) as the symmetry-axis. For particles larger 
than 4 micron in diameter the ion drag can be ignored in the 
vertical force balance lH3l . so that we have 



- m D g - q D 



d<j>(r, z) 



dz 



F th = 0, 



(2) 



where is the dust mass, which we assume to know, g 
the magnitude of the gravitational acceleration, qp the dust 
charge. The induced thermophoretic force is given by 



Fth 



32 n T R 2 
15 vt 



VT, 



gas i 



(3) 



where kt = 0.0177 is the argon heat conduction coefficient, R 
the particle radius, Vt the thermal velocity of the neutral atoms 
and VTg as the gradient in temperature of the gas between the 
electrodes. For the radial confinement we have 



qD- 



d< j>{r, z) 
dr 



= —kr. 



(4) 



By using partial differentiation and combining these equations 
we find the equation for k, 



k = 2c [m D g - F th ] 



(5) 



which thus depends on the steepness of the potential well and 
the applied electrode temperature. 

Assuming that our crystal is single-layer and consists of TV 
particles with inter-particle spacing A (which depends on r) 
interacting through an inter-particle potential V(A), it can be 
shown [13] that the pressure in the crystal layer can be written 
as 



P = 



1 d[3NV(A)} _ V3 dV(A) 
'N d(VSA 2 /2) ~ ~A dA 



(6) 



Using appropriate boundary conditions, the theoretical edge of 
the crystal can be shown to be 



Ao 
Xd 



(7) 



where Ao is the central inter-particle spacing. The actual 
observable outer edge of the crystal, Rm, is given by Rm = 
Roo — Am\/3/2, with Am the inter-particle spacing at the 
edge of the crystal. The total number of particles in the layer, 
N can be derived as 



N = "TIT — T" + T-) V ( A o)- 

kA \A X D J v uy 



(8) 



IEEE TRANSACTIONS ON PLASMA SCIENCE, VOL. XX, NO. XX, XXXXX 2010 



3 



Combining the above two equations and solving for V(Ao), 
we can obtain an equation for the Debye length, 



Xd = A 



A-S 
3S-A 



(9) 



where we have defined the total crystal surface area as 
A = ttR 2 ^ and the total surface area covered by N Wigner- 
Seitz cells measured at the center to be S = V3AqN/2 
After determining the Debye length, the dust charge can be 
calculated from 



Qd 



47re AoA;i?2 ) 



\3(3 + fE#)exp(^f)" 



(10) 



Once we have independently determined fco, the constant of 
the radial confining force without applied thermophoresis, we 
can calculate c and then we can calculate k for any additionally 
applied thermophoresis from equation (0). After measuring 
the observables Rm, N, Ao, and Am, we can determine 
the Debye length and the dust charge from equation © and 
equation (TTOb , respectively. When we obtain the levitation 
height of the crystal from a side-view picture, we can then use 
the force balance to reconstruct the electric field as a function 
of height. 

A. Measuring k without thermophoresis, ko 

The common assumption concerning the radial confinement 
resulting from cutouts or other adaptations to the lower elec- 
trode is that it can be described by a potential well varying 
quadratically with the distance from the center. Even though 
this has been shown for specific experiments fl6lh this is not 
a priori clear for our case. 
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Fig. 2. The natural logarithm of the distance from the position in the final 
frame plotted versus time for four selected particle trajectories. The dashed 
lines are fits (R 2 > 0.99) according to the solution of equation (12\ . 

In order to determine the radial restoring force, we ran the 
discharge at high pressure (1 Torr in this case) in order to 
minimize any oscillatory motion of the dust particles, without 
heating or cooling of the lower elecotrode, so that F t h = 0, 
and dropped the dust particles just outside of the cutout. 
We then monitored the particles while they moved inwards 
into the potential well of the cutout. Assuming that in this 
case the forces acting on the particles are the electrostatic 



force (determined by the potential well) and the gas drag, the 
equation of motion for these particles can be written as: 



tudx + ax + kox = 0, 



(11) 



with a the Epstein neutral drag friction coefficient ifTTl . We 
assume that the particles are monodisperse, so that they all 
have the same mass, mo- 

The general solution to equation (ITTI) for the overdamped 
case (when a 2 > 4m£>fco) is given by a sum of two exponen- 
tials, 



x(t) = Aexp(Cif) + £exp(C 2 t), 
with Ci 5 2 gi yen by 



C, 2 = - -—±\ (—) -4*L 
2 rn D y \m D J m D 

Taking their sum and difference we find 



Cx+C 2 = 

m D 



(12) 



(13) 



C\ — C2 



(14) 
(15) 



Hence, fitting trajectories of particles moving in the potential 
well provides us with C\ and C2, from which we derive a/rriD 
and ko/rriD- 

Figure [2 shows four selected particle trajectories. Plotted is 
the natural logarithm of the distance from the final position in 
the last frame. The dashed lines are fits according to equation 
(TTZl) . The fits are very reliable (R 2 > 0.99). From the fits, we 
determine a/mo = 14±1.9s _1 and ko/mo = 6.4±0.87 s -2 . 
For our 8.9 micron MF particles, we have mrj = 5.55 x 10 -13 
kg, so that we have a = 7.8 ± 1.1 x 10 -12 kg/s, which is in 
reasonable agreement with results reported in lfl8l . For fco we 
thus find k = 3.6 ± 0.5 x 10~ 12 kg/s 2 . 

IV. Determination of the Debye length, charge, 

AND ELECTRIC FIELD 

Now that we have determined fco f° r me 1-inch diameter 
cutout in our modified GEC cell, we can take pictures of 
the dust crystals formed with different applied temperature 
gradients. From side-view pictures the levitation height was 
determined, and from top-view pictures the number of par- 
ticles, crystal radius, and central inter-particle distance were 
determined using graphics software lfT9l . 

Figure [3] shows a top-view image of a typical crystal 
levitated in the discharge. This crystal was formed while 
the lower electrode was cooled to 0° C. It contains roughly 
N=1000 particles, has a radius of Roo= 8.2 mm, and the central 
inter-particle distance was measured to be A = 386 fim. This 
yields A = 2.4 x 10~ 4 m 2 and S = 1.3 x 10~ 4 m 2 , so that 
from equation ©, we obtain A^> w 390/im. 

Since the thermophoretic force acts in the same direction as 
gravity when the electrode is cooled, k increases with decreas- 
ing lower electrode temperature. For the applied temperature 



IEEE TRANSACTIONS ON PLASMA SCIENCE, VOL. XX, NO. XX, XXXXX 2010 



4 




Fig. 3. A top- view of a dust crystal composed of about 1000 particles 8.9 
fim in diameter, at an electrode temperature of 0° C, applied power of 2 W, 
and a pressure of 200 mTorr. The crystal has a radius of 8.2 mm and a central 
inter-particle distance of 386 /im. The image has been processed to improve 
the visibility of the dust particles. 



gradient of 906 K/m, we find that k (0° C) is 5.3 x 10" 12 kg/s 2 . 
Using this and the value for in equation (ITOb , we find a dust 
charge of qp ~ -1.3 xlO 4 e. Using the vertical force balance, 
this results in an electric field of E = (m^ — F t h)/qD ~ 
4600 V/m. From side-view pictures, like the one shown in 
figure |U the levitation height of the upper crystal-layer was 
determined to be 5.8 mm. 
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Fig. 5. The levitation height of the dust crystal versus temperature of the 
lower electrode, determined from side-view images. The dashed line is an 
exponential fit to the data. 
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Fig. 6. The Debye length, dust charge, and electric field plotted versus 
height above the electrode, determined from top- view and side-view pictures. 
The electrode temperature was varied between 0° C and 60° in steps of 10 
degrees. 



Fig. 4. A side-view of a dust crystal layer composed of 8.9 /mi diameter 
particles. It can be seen that the crystal levitated in the sheath is not completely 
single-layer. The gray area in the upper-right corner is glow from the plasma. 
The figure has been processed to improve the visibility of the dust particles. 



Repeating these calculations for different applied lower 
electrode temperatures, the Debye length, dust charge, and 
electric field were obtained as a function of the temperature. A 
plot of the dust levitation height versus applied temperature is 
shown in figure The fit shows an exponential dependence of 
the levitation height with the temperature. The Debye length, 
dust charge, and electric field are shown in figure 

As shown, the Debye length varies between 200 and 360 
/im, having its smallest values higher in the plasma. This 
seems reasonable, since the plasma density is expected to 
increase farther away from the electrode, and ex y/T/n. 
Since the ion temperature is closely coupled to the neutral gas 
temperature, it has to be mentioned that when the screening 
is mainly determined by the ions, that the increase in gas 
temperature could lead to an increase in ion Debye length. 
The data suggest that the decrease in plasma density with 
height is more important. The measured values for the Debye 
length are in agreement with values reported in the same GEC 
cell (roughly 250 /im for a discharge at 200 mTorr and 5 W 



power), measured using a different technique Q, and also with 
the observed inter-particle distances, which in these crystals 
typically are a few times the Debye length |4|. 

The dust charge becomes continuously less negative from 
approximately -1.3xl0 4 e to -5.5xl0 3 e with increasing height 
above the lower electrode. These values are lower than the 
value predicted by OML theory [20], which is in agreement 
with previous results, for instance those reported in EH . 

The electric field decreases approximately linearly for 
heights above 6.5 mm. The obtained vertical electric field in 
the sheath varies from -4600 V/m (closer to the electrode) 
to -1380 V/m (farther up in the plasma). The drift velocity 
Vd = H+E for argon ions (with mobility of roughly 1 m 2 /Vs 
ll22l ) then varies from 4600 m/s to 1380 m/s. According to the 
Bohm criterion, the ions are entering the sheath at the Bohm 
velocity, vb = \/kBT e /m+ m/s, which, for typical electron 
temperature of 3-5 eV |14], corresponds to 2700-3500 m/s. 
We conclude that we are probing the plasma both above the 
Bohm point, consistent with fT3lL as well as just below the 
Bohm point. 
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V. Discussion and conclusions 

There are three main assumptions in the analysis of our 
measurements. First, we assume that the particles are radially 
confined in a parabolic potential well provided by the cutout 
in the electrode, second, that the ion drag force does not play 
a role in the vertical force balance for the 8.9 micron diameter 
particles levitated in the sheath, and third that the equation of 
state derived for a single-layer dust crystal and the resulting 
equations for the observables are valid in our experiment. 

The trajectories of particles moving into the potential well 
seem to behave very much like an overdamped harmonic 
oscillator with a restoring force in Hooke's form, so that 
F = —dU(r)/dr — —kr, hence U(r) ex cr 2 . There was 
no indication of higher order terms playing a role in the 
trajectories. This might not be true much closer to the cutout 
edge, but plays no role for our current analysis. 

Analysis in fl3l has shown that the ion drag force needs 
only to be taken into account for particles smaller than 4 fim 
in diameter. Since then, however, the understanding of the ion 
drag has improved, for instance where non-linear scattering, 
ion-neutral collisions, and anisotropic screening due to ion 
flow are concerned, which might be important in the sheath 
l23l - l25ll . We cannot completely rule out the role of the ion 
drag, but the values obtained for the charge, Debye length and 
the electric field are consistent with values obtained from the 
literature, so that we expect the role of the ion drag to be small 
in our case. 

Side-view pictures show that the observed crystals were 
not completely single-layer, but that layer splitting occurs, 
especially at higher electrode temperatures. Defects induced 
by this process are a source of free energy in the crystal and 
therefore the equation of state might not be correct. However, 
these defects mainly affect the crystal locally. The number of 
particles and the crystal radius are global observables, which 
are probably not affected by these defects. The compression 
of the crystal in the parabolic potential well was still observed. 
When measuring inter-particle distances we took averages over 
many different inter-particle distances at various locations. All- 
in-all, we think that the local deviations from the pure single- 
layer crystal state have not affected our results strongly. 

The systematic error due to the assumptions in our model 
are believed to be small, but a word about the precision in 
obtaining our observables is then required. Surprisingly, the 
observable hardest to obtain was the total number of particles, 
both due to the fact that the crystals did not always fit in 
our field-of-view while maintaining the resolution necessary 
to resolve the dust grains, as well as due to the layer- splitting. 
We estimate the error to be roughly 10-15%. The distances 
measured from the top- view pictures were much more precise. 
We estimate that any measurement of the radius is off at most 
by 3-5 pixels, given our method of fitting a circle to the crystal. 
This corresponds to roughly 50-80 micron, which for an 8 
mm radius crystal is an error of at most 1%. The inter-particle 
distances could be off by a similar number of pixels, but the 
results are the averages of many measurements. The error is 
anticipated to not exceed 5% because of this. We see that the 
largest error in our experiment is actually presented by the 



determination of fco, which is roughly 14%. 

An interesting observation, which we could not find in the 
literature, was that a measurable change in the natural DC bias 
on the lower electrode occured. This bias became less negative 
for electrode temperatures above room-temperature, changing 
from -13 V when the electrode was at room-temperature to 
roughly -10.5 V at +60° C. When the electrode was cooled 
below room-temperature, we did not observe the opposite 
effect as clearly. Since a decrease in the negative DC bias 
causes dust particles to levitate closer to the lower electrode, 
negating the effect of the applied temperature, we maintained 
the DC bias at -13 V for all temperatures using an external 
power supply. This insured that the DC bias did not play a 
role in the levitation of the dust. 

In conclusion, we have succesfully used an entire dust 
crystal layer as a probe for the electric field in the sheath 
in front of the powered electrode in a modified GEC cell, 
by employing thermophoresis. We obtained the Debye length, 
the dust charge and the electric field by simply measuring 
observables from top-view pictures and by measuring the lev- 
itation height from side-view pictures. This method provides a 
single, non-perturbative and practical method, which does not 
require complex alterations to dusty plasma experiments. We 
have shown that 8.9 fim diameter MF particles are mostly 
suspended on the plasma side of the Bohm point in our 
discharge at typical discharge settings. 
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